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We examine how the new forthcoming Earth gravity models from the CHAMP and, 
especially, GRACE missions could improve the measurement of the general relativis- 
tic Lense— Thirring effect according to the various kinds of observables which could be 
adopted. In a very preliminary way, we use the recently released EIGEN2 CHAMP-only 
and GRACEOIS GRACE— only Earth gravity models in order to assess the impact of the 
mismodelling in the even zonal harmonic coefficients of geopotential which represents 
one of the major sources of systematic errors in this kind of measurement. 



1. Introduction 

The post-Newtonian gravitomagnetic Lense-Thirring efFect^'^^'^, in the slow motion 
and weak field approximation of the General Theory of Relativity, consists of secular 
precessions of the longitude of the ascending node fl and the argument of pericentre** 
uj of the orbit of a test particle freely falling in the gravitational field of a central 
spinning object of mass M and proper angular momentum J 

f2LT = 3 , (1) 

c^a'^(l — e ) 2 

QGJcosi , . 

^^'^ = — ttt; — iT^' 

where a, e and i are the semimajor axis, the eccentricity and the inclination, re- 
spectively, of the orbit of the test particle and G is the Newtonian gravitational 
constant; c is the speed of light in vacuum. 



1.1. The LAGEOS LAGEOS II Lense Thirring experiment 

The first attempts to measure such effect in the gravitational field of Earth are 
due to Ciufolini and coworkers which analysed the data of the orbits of the exist- 
ing LAGEOS and LAGEOS II geodetic satellites collected with the Satellite Laser 

^In the original paper by Lense and Thirring the longitude of the pericentre to = £7 + oj is used 
instead of u). 
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Ranging (SLR) technique*. The observable used is the following combination^ of 
the orbital residuals of the nodes of LAGEOS and LAGEOS II and the perigee of 
LAGEOS II 

^fjLAGEOS ^ ^^sh'^AGKOS II ^ ^2^^^^^°^ " ~ 60.2mlT, (3) 

with 

ci = 0.304, ci = -0.350, (4) 
The quantity /^lt is a solved-for least square parameter which is in Newtonian 

Table 1. Orbital parameters of LAGEOS, LAGEOS II, LARES, SI and S2. 



Orbital parameter 


LAGEOS 


LAGEOS II 


LARES 


SI 


S2 


a (km) 


12270 


12163 


12270 


12000 


12000 


e 


0.0045 


0.014 


0.04 


0.05 


0.05 


i (deg) 


110 


52.65 


70 


63.4 


116.6 



mechanics and 1 in the General Theory of Relativity. The gravitomagnetic signature 
is a linear trend with a slope of 60.2 milliarcseconds per year (mas yr~^ in the 
following). The claimed total accuracy is of the order of 20% 30%. In Table 1 we 
quote the orbital parameters of the existing or proposed LAGEOS-like satellites 
useful for general relativistic tests. 

The main sources of systematic errors in this experiment are 

• the unavoidable aliasing effect due to the mismodcUing in the classical sec- 
ular precessions induced on and u) by the even {i = 2, 4, 6,...) zonal 
(to = 0) coefficients^ of the multipolar expansion'^ of the terrestrial 
gravitational field 

• the non-gravitational perturbations affecting especially the perigee of LA- 
GEOS II. Their impact on the proposed measurement is difficult to be 
reliably assessed 

It turns out that the mismodelled classical precessions due to the first two even 
zonal harmonics of the geopotential J2 and J4 are the most insidious source of error 

in the measurement of the Lense-Thirring effect with the nodes of LAGEOS and 
LAGEOS II and the perigee of LAGEOS II only. It is so also because of their 
secular variations J2, J4, which can be accounted for by an effective J2 coefficient 
jf ~ j2-F0.371 j4-f 0.079j6 + 0.006j8-0.003jio..., and of the fact that the Lcnsc- 
Thirring effect is embedded in their recovered values released by any Earth gravity 



"Jl = —CiQ, ( = 2, 4, 6, .... The unnormalized Stokes coefficients of degree i and order m 

can be obtained from the normalized Stokes coefficients according to C^^n = ^tmCtm where 

^ _ \ (.2e+l){2-5a^)(l-m)\ 

"^It accounts for the oblateness of Earth generated by its diurnal rotation. 
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models since the gravitomagnetic acceleration is a priori included in the force models 
used in obtaining, among other things, the even zonal coefficients. The combination 
of eq.|j2Jl is insensitive just to J2 and J4. According to the full covariance matrix of 
the EGM96 gravity model^, the error due to the remaining uncancelled even zonal 
harmonics amounts to almost^ 13%. However, if the correlations among the even 
zonal harmonic coefficients are neglected and the variance matrix, in a root-sum- 
square fashion, is used*^, the error due to the even zonal harmonics of geopotential 
amounts to^ 45%. With this estimate the total error of the LAGEOS-LAGEOS 
II Lense-Thirring experiment would be of the order of 50%. If the sum of the 
absolute values of the individual errors is considered, an upper bound of 83% on 
the systematic error due to the even zonal harmonics is obtained; the total error 
in the LAGEOS-LAGEOS II Lense-Thirring experiment would, then, amount to 
almost 100%. It must be considered that, since in EGM96 a strong correlation 
among the various recovered even zonal coefficients does occur, the approach of the 
sum of the absolute values of the individual errors would, probably, yield a more 
realistic estimate of this kind of systematic error. 

1.2. The supplementary orbital planes option 

The originally proposed LARES mission^ consists of the launch of a LAGEOS-type 
satellite-the LARES-with the same orbit of LAGEOS except for the inclination i of 
its orbit, which should be supplementary to that of LAGEOS, and the eccentricity 
e, which should be one order of magnitude larger. The choice of the particular value 
of the inclination for LARES is motivated by the fact that in this way, by using as 
observable the sum of the nodes of LAGEOS and LARES 

E5f2 = f2LAGEOS^^j=^LARES^ (5-, 

it should be possible to cancel out, to a high degree of accuracy, all the contributions 
of the even zonal harmonics of the geopotential, which depends on cosi, and add up 
the Lense-Thirring precessions which, instead, are independent of i. The systematic 
error due to the even zonal coefRcients of geopotential would range from^° 0.8% to 
4.6% for an uncertainty in ilares up to 1 deg with respect to its nominal values 
of Table 1, according to the variance matrix of EGM96 used in a root-sum-square 
fashion. If, instead, the sum of the absolute values of the individual errors is consid- 
ered, the systematic error due to the even zonal coefficients of geopotential would 
amount to almost 2% for jlares = 70 deg and to 10% for a departure of ilares 



"^Such approach is considered more realistic by some authors* because nothing would assure that 
the correlations among the even zonal harmonics of the covariance matrix of the EGM96 model, 
which has been obtained during a multidecadal time span, would be the same during an arbitrary 
past or future time span of a few years as that used in the LAGEOS-LAGEOS II Lense-Thirring 
experiment or in the proposed LAGEOS-LARES mission. Moreover, in the same paper it is 
claimed that a correct and reliable evaluation of the impact of the non-gravitational perturbations 
of the perigee of LAGEOS II on the proposed measurement is rather troublesome. 
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of 1 deg from its nominal value. The impact of the non-gravitational perturbations 
would be far less important because the nodes of the LAGEOS-like satellites are 
almost insensitive to them. 

Since the Lense-Thirring precession of the pericentre depends on cos i and the 
classical precessions due to the even zonal harmonics of geopotential depend on 
cos^ i and odd powers of sinz, also the difference of the perigees^^'^^ 

Auj = Slu^'^ - (6) 

of a pair of satellites with identical orbits in supplementary orbital planes could 
be used, in priciple, in order to measure the Lense-Thirring effect. It has been 
shown that, for a pair of LAGEOS-like satellites SI and S2 (see Table 1) with the 
same physical properties of the LAGEOS satellites, the total error would not be 
less than^^ 5%, mainly due to the non-gravitational perturbations. 

1.3. Other approaches 

In order to cope with the problem of the impact of the unavoidable orbital injection 
errors in the inclination of LARES on the systematic error due to the even zonal 
harmonics of geopotential, the following observable has also been proposed^'^ 

^j^LAGEOS ^ CK^iil^AGEOS II ^ ^^Sn^^ARES ^ ^g^^LARES _ gl.g^LT, (7) 

with 

Cl ~ 3 X 10"^ C2 - 9.9 X lO-\ C3 - 1 X lO^'l (8) 

It would allow to neglect the problem of the correlations between the various even 
zonal coefficients of geopotential. Moreover, the secular variations in the even zonal 
harmonic coefficients of geopotential would not affect such observable. The same 
would also hold for the bias of the Lense-Thirring effect itself on the even zonal 
harmonics because it is mainly concentrated in J2 and J4. Notice that such an 
observable does not include the perigee of LAGEOS II. 

Another possible way to overcome the uncertainties related to the error budget 
evaluation due to the non-gravitational perturbations which affects the perigee of 
LAGEOS II by using the existing satellites is the following combination^^ 

<5liLAGEOS ^ ^^^j^jLAGEOS II _ 48.2^lT, Ci = 0.546 (9) 

A similar approach is implicitly presented in ref.^^. The observable of eq.® would 
be affected neither by the secular variations in the Jg coefficients nor, at least to a 
certain extent, by the a priori Lense-Thirring bias of the even zonal harmonics. 

1.4. The role of the various Earth gravity models 

Here we wish to discuss the impact of the new, forthcoming Earth gravity models 
from the CHAMP^^ and GRACE^^ missions, which should greatly improve our 
knowledge of the terrestrial gravitational field, on the different combinations used 
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or proposed for the measurement of the Lense-Thirring effect. Very recently, the 
EIGEN2i^ GRACEOIS^^ and GGM01C/S'= models from the CHAMP and GRACE 
only data, respectively, have been released. Although very preliminary and lacking 
of reliable validation and calibration tests, they can be used in order to get an idea 
of which the improvements in the Lense-Thirring tests could be. In Table 2 we 
quote the sigmas of the normalized even zonal Stokes coefficients cr^^^ for EGM96, 
EIGEN2 and GRACEOIS. Note that for the GRACE models the sigmas are not the 
formal errors but are approximately calibrated. For EIGEN2 the formal errors are 
released: it is likely that they are rather optimistic, at least for the low degree even 
zonal harmonics up to £ = 20^^ . It should be noted that J2 is known much less 

Table 2. Errors in the even zonal normalized Stokes coefficients f^c^^ for various Earth gravity 
models up to degree £ = 20. The GGMOIC model is very similar to GRACEOIS, apart from 
(Tg^^ = 1.1677 X 10~" and a^^^ = 1.7874 X 10"". 



i 


EGM96 




EIGEN2 


GRACEOIS 


2 


3.561x10- 


11 


1.939x10" 


ii 


3.106x10" 


-10 


4 


1.042x10" 


■10 


2.230x10" 


11 


2.713x10" 


'11 


6 


1.449x10" 


■10 


3.136x10" 


11 


9.568x10" 


-12 


8 


2.266x10" 


■10 


4.266x10" 


11 


6.428x10" 


-12 


10 


3.089x10" 


10 


5.679x10" 


11 


7.180x10" 


-12 


12 


4.358x10" 


10 


7.421x10" 


11 


5.022x10" 


-12 


14 


5.459x10" 


10 


9.561x10" 


11 


4.376x10" 


.12 


16 


5.313x10" 


10 


1.216x10" 


10 


3.805x10" 


-12 


18 


4.678x10" 


10 


1.532x10" 


10 


4.412x10" 


-12 


20 


4.690x10" 


■10 


1.914x10" 


10 


3.571x10" 


-12 



accurately in GRACEOIS than in EGM96 and EIGEN2. 

2. Combinations with the existing or proposed geodetic satellites 
of LAGEOS-type 

The classical secular precessions of the node and the perigee, for a given even zonal 
harmonic of degree £, are proportional to R a ^ ^ > , where R is the Earth mean 
equatorial radius. Then, the relevant classical secular precessions of the node and 
the perigee of the LAGEOS satellites, due to their high altitude (see Table 1), are 
those induced just by the first few five-six even zonal harmonics of geopotential 
in the sense that the error in the Lense-Thirring measurement due to geopotential 
does not change if the even zonal harmonics of degree higher than £ = 10-12 are 
neglected in the calculations. 



'^It can be retrieved on the WEB at |http : / / www. csr . utexas . edu/ grace/gr avityT] The GRACE- 
only GGMOIS model was combined with the TEG-4 information equations (created from historical 
multi— satellite tracking data; surface gravity data and altimetric sea surface heights) to produce 
the preliminary gravity model GGMOIC. 
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2.1. The full range even zonal harmonics observables 

Let us recall that the sum of the nodes Eli of LAGEOS and LARES (and of SI 
and S2), and the difference of the perigees Aoj of SI and S2 do not cancel any even 
zonal harmonic coefficients of geopotential due to the unavoidable orbital injection 
errors; in the case of the originally proposed LAGEOS-LARES mission this non 
perfect cancellation of the effects of the even zonal harmonics would occur even 
if the orbital parameters of LARES were exactly the same as in Table 1 because 
Clares = 0.04, while clageos — 0.0045. This impHes that it is of the utmost 
importance that the new gravity models yield notable improvements especially for 
J2, J4 and, to a lower extent, Jq and Jg which would affect such observables and 
which are the most relevant sources of systematic errors for the Lense-Thirring 
precessions of the nodes and the perigees of each satellite. This is also the case 
for the combination of eq.® which involves the node of LAGEOS and the node of 
LAGEOS II; it cancels out the first even zonal harmonic J2 and is affected by the 
remaining J4, Jg,... 

Let us consider the sum of the nodes of LAGEOS and LARES. By assuming an 
orbital injection error in zlares up to 1 deg, it turns out that the variance matrix 
of EIGEN2 up to degree £ = 20 yields an error due to the even zonal harmonics of 
geopotential spanning from 0.2% to 1.6%. GRACEOIS, instead, induces an error 
which ranges from 0.2% to 24%. The range error induced by GGMOIC is 0.2%-0.9%. 
Let us recall that the error due to EGM96 ranges from 2% to*' 10%. These numbers 
can be explained if we consider that for large departures from the condition of ideal 
supplementary orbital planes the contributions of all the even zonal harmonics-and, 
especially, those of lower degrees-affect the systematic error due to geopotential. 
From Table 2 it can be noticed that J2 is far better known in EIGEN2 than in 
GRACEOIS. An analysis of the contributions of the various even zonal harmonics in 
the building up of the systematic error due to geopotential shows that for departures 
of Glares of 1 deg from its nominal value of Table 1 the most important role is played 
just by J2. This also explains the better results obtainable with GGMOIC. 

The combination of eq.© is affected by the uncancelled even zonal harmonics J4, 
Jg,... at a 177%, 22%, 21% and 14% level according to the sum of the absolute values 
of the errors of EGM96 and to the variance matrices of EIGEN2, GRACEOIS and 
GGMOIC, respectively. The small difference between the EIGEN2 and GRACEOIS 
results can be explained by noticing that J2 does not affect eq. ^ , the error in J4 is 
shghtly better in EIGEN2 than in GRACEOIS and the accuracy of the GRACEOIS 
determinations of Je, Js,.-- is only slightly better than that of EIGEN2. The better 
results for GGMOIC, especially for J4, explain the result obtained with such Earth 
gravity model. 



'Note that for EIGEN2, GRACEOIS and GGMOIC a root-sum-square calculation with their 
variance matrices yields realistic estimates because, contrary to EGM96, their covariance matrices 
are almost diagonal. 
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2.2. The partial range even zonal harmonics observables 

If, instead, we consider the suitably designed combinations of N orbital residuals 
of the LAGEOS satellites which cancel out the first N — 1 even zonal harmonics^ 
the situation is different in the sense that the improvements of the remaining higher 
degree even zonal harmonics become relevant in order to reduce the systematic error 
due to the even zonal harmonics of geopotential which, for > 3 is well below 1% 
also with EGM96 (RSS calculation). 

Let us consider the currently used observable of eq. Q which cancels out J2 and 
J4. According to the variance matrices of EGM96, EIGEN2 and GRACEOIS up 
to degree £ = 20 the systematic error due to geopotential amounts to 83% (sum 
of the absolute values of the errors), 9% and 2% (root-sum-square calculations), 
respectively. This shows the importance of reducing the uncertainties in Jg and, to 
a lower extent, Jg. 

The combination of eq.ljZj) is affected by the remaining even zonal harmonics Jg, 
Jio,... far less than 1%. Indeed, for a maximum orbital injection error in ilares of 
1 deg the systematic error due to geopotential amounts to 2.2%-3% according to 
the sum of the absolute values of the errors in the variance matrix of EGM96 up to 
degree £ = 20, to 0.08%-0.1% according to the variance matrix of EIGEN2, used in 
a root-sum-square fashion, up to degree £ = 20 and to 0.007%-0.01% according to 
the variance matrix of GRACEOIS, used in a root-sum-square fashion, up to degree 
i = 20. This can be explained by noticing from Table 2 that the higher degree even 
zonal harmonics are far better determined in GRACEOIS (and GGMOIC/S) than 
in the other Earth gravity models. 

3. Combinations with the other existing geodetic satellites 

In principle, a very appealing possibility would be the use of combinations which 
include also the orbital elements of the other existing geodetic satellites like Ajisai, 
Starlette, Stella, WESTPACl, ETALONl, ETAL0N2, apart from those of LA- 
GEOS and LAGEOS II; this option has been studied in some previous works^^. In 
Table O the orbital elements of the other existing SLR satellites are reported. For 

Table 3. Orbital parameters of Ajisai, Starlette, Stella, WESTPACl, ETALONl and ETALON2 



Orbital parameter 


Ajisai 


Starlette 


Stella 


WESTPACl 


ETALON 1 


ETALON2 


a (km) 


7870 


7331 


7193 


7213 


25498 


25498 


e 


0.001 


0.0204 








0.00061 


0.00066 


i (deg) 


50 


49.8 


98.6 


98 


64.9 


65.5 



example, one might think about some combinations including only the nodes of the 
LAGEOS satellites and of some of the other geodetic SLR satellites. Unfortunately, 



^We assume N > 2. 
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it turns out that this way is unpracticable because the other satelUtes orbit at lower 
altitudes than LAGEOS and LAGEOS II (see Table OJ, so that they are sensitive 
to much more even zonal harmonics, especially Starlette, Stella and WESTPACl. 
Also EIGEN2 and GRACEOIS do not alter this situation in a significative way. In- 
deed, it turns out that, while the systematic error due to the even zonal harmonics 
of geopotential of a combination which includes the nodes of LAGEOS and LA- 
GEOS II, the perigee of LAGEOS II and the node of Ajisai is sensitive only to the 
harmonics up to degree £ = 14, according to the variance matrix of GRACEOIS, for 
another combination similar to the previous one, with the node of Starlette instead 
of the perigee of LAGEOS II, the error induced by geopotential is sensitive well 
up to degrees £ = 40-50. Moreover, the calculation of the classical secular nodal 
precessions for degrees as high as £ = 40 is rather unreliable because of instabilities 
of the numerical results. However, it may be interesting to consider the following 
combination 

with ci = 4.174, C2 = -2.705, C3 = 1.508, C4 = -0.048. According to a Root- 
Sum-Square calculation with the variance matrix of GGMOIC up to^ £ = 42 the 
impact of the remaining even zonal harmonics of degree £ > 10 amounts to 21.6 mas 
yr^^ which yields a 37.6% percent error due to geopotential in the measurement 
of the Lense-Thirring effect with eq. fllHl . The upper bound due to the sum of the 
absolute values of the individual errors amounts to 123%. If the future GRACE- 
based gravity solutions will improve the high degree (Jio, J12, Ju, •■•) even zonal 
harmonics more than the low degree (J2, J4, Jq, Js) ones, the combination of 
eq. (|10|l might deserve some interest in alternative to that of eq. . 

In conclusion, also with the new forthcoming Earth gravity models, it is unlikely 
that the combinations involving the orbital elements of the other existing geodetic 
satellites as well will represent a real alternative to those which use only the existing 
LAGEOS satellites and the proposed LARES. If, instead, we refer to the existing 
LAGEOS satellites only, it may happen that the inclusion of the nodes of the other 
geodetic satellites could become a viable possibility. 

4. Conclusions 

The choice of the optimal observable for measuring the Lense-Thirring effect is 
dictated by the need of getting a compromise among the reduction of the system- 
atic gravitational error due to geopotential, the systematic non-gravitational error 
induced by the perturbations affecting especially the perigees of the geodetic SLR 
satellites and the reliability of the SLR data of the satellites orbits. 



''It has been checked that the error due to the even zonal harmonics remains stable if other even 
zonal harmonics are added to the calculation. Moreover, it turns also out that, up to ^ = 42, there 
are no appreciable fluctuations in the calculated classical secular precessions. 
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In regard to those requirements, it turns out that even with the new, forthcoming 
Earth gravity models the multisatelhte approach involving the orbital elements of 
the existing SLR satellites other than the LAGEOS ones should not be competitive 
with the observables built up with the LAGEOS-type satellites only. 

Among them, notice that the total accuracy obtainable with the originally pro- 
posed combination by Ciufolini, which involves the nodes of both LAGEOS and LA- 
GEOS II and the perigee of LAGEOS II, strongly depends on the correct evaluation 
of the impact of the non-gravitational perturbations on the perigee of LAGEOS II 
because it is arguable that the error due to geopotential could be reduced to 1% 
or, perhaps, less in the near future. On the contrary, for the combination proposed 
explicitly in ref.^^, which involves the nodes of LAGEOS and LAGEOS II, the im- 
provements in the Earth gravity field even zonal coefficients, especially J4, Jq, Jg, 
will be of crucial importance in reducing the total error to less than 10% because 
the impact of the non-gravitational perturbations on the error budget is quite neg- 
ligible. Also a combination involving the nodes of LAGEOS, LAGEOS II, Ajisai, 
Starlcttc and Stella might become interesting, provided that notable improvements 
in the higher degree even zonal harmonics Jio, J12, ... will occur. 

For the observables involving the proposed LARES, the multisatelhte combina- 
tion including also the elements of LAGEOS and LAGEOS II would not require 
any particular improvements in the knowledge of the higher degree even zonal har- 
monics of geopotential because the systematic error induced by them, according to 
the present-day level of knowledge of the terrestrial gravitational field, would be 
sufficiently small. 
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